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ABSTRACT 

We study the abundance of satellite galaxies as a function of primary stellar mass us- 
ing the SDSS/DR7 spectroscopic catalogue. In contrast with previous studies, which fo- 
cussed mainly on bright primaries, our central galaxies span a wide range of stellar mass, 
10 7 ' 5 ?C M^ 1 /Mq ^ 10 11 , from dwarfs to central cluster galaxies. Our analysis confirms 
that the average number of satellites around bright primaries, when expressed in terms of 
satellite-to-primary stellar mass ratio (m^'/Mf 11 ), is a strong function of M* rl . On the 
other hand, satellite abundance is largely independent of primary mass for dwarf primaries 
(M* rl < 10 10 Mq). These results are consistent with galaxy formation models in the ACDM 
scenario. We find excellent agreement between SDSS data and semi-analytic mock galaxy 
catalogues constructed from the Millennium-II Simulation. Satellite galaxies trace dark mat- 
ter substructure in ACDM, so satellite abundance reflects the dependence on halo mass, 1/200, 
of both substructure and galaxy stellar mass (M*). Since dark matter substructure is almost 
scale-free, the dependence of satellite abundance on primary mass results solely from the well- 
defined characteristic mass in the galaxy mass-halo mass relation. On dwarf galaxy scales, 
where models predict a power-law scaling, M» oc M^qq, similarity is preserved and satel- 
lite abundance is independent of primary mass. For primaries brighter than the characteristic 
mass of the M*-M2oo relation, satellite abundance increases strongly with primary mass. Our 
results provide strong support for the steep, approximately power-law dependence of dwarf 
galaxy mass on halo mass envisioned in ACDM galaxy formation models. 

Key words: galaxies: dwarf - galaxies: formation - galaxies: haloes 
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1 INTRODUCTION 

Matching the galaxy luminosity function in the ACDM sce- 
nario requires that the stellar mass of galaxies, M», should 
vary strongly with the viria{JJ mass, A/200, of their surrounding 
dark matter halos. This exercise implies that the "efficiency" of 
galaxy formation, as measured by the ratio M*/M2oo, decreases 
steadily toward both smaller and larger masses from a maxi- 
mum at M200 ~ , 1 ,° 12A f ' 7 : <fYanget al.ll2003l ; IVale & OstrikeJ 
, 2004 IShankar et al.l l2006UZheng et ai.M2007l; iKravtsoyl |2qTcF 
Moster et alj|20ld : IBehroozi et alJl2oioF IGuo et alj I2010L 1201 lh . 
On the scale of dwarf galaxies (M* < 10 10 Mq) these models 
require a near power-law dependence, M» oc Mfpo , in order to re- 
produce observations of faint objects (however, see IBehroozi et al.l 



Virial quantities are defined at the radius from the center of each halo 
where the mean enclosed density equals 200 times the critical density of 
the Universe and are identified by a 200 subscript. Units assume a Hubble 
constant of Hq = 73 km s _1 Mpc -1 unless otherwise specified. 



|2012|) . Such a steep M* -A/200 relation implies that dwarfs differ- 
ing by as much as three decades in stellar mass (e.g., from the For- 
nax dwarf spheroidal to the Large Magellanic Cloud) inhabit halos 
spanning just over one decade in virial mass. Furthermore, very few 
galaxies exceeding 10 6 Mq are expected to have halos with virial 
masses below 10 10 Mq jFerrero et a"i]|201 ll) . 

These predictions have been recently challenged by a series of 
observations, including (i) the lack of a characteristic velocity at the 
faint-mass en d of blind HI surveys (expected if most dwa rfs live in 
similar halos, IZwaan et alJl201Cc IPapastergis et alj|201ll) : and (ii) 
the low virial mass (substantially below 10 10 Mq) inferred from 
dyna mical data for the dwarf sph eroidal companions of the Milky 
Way (B oylan-Kol phin et al]|2012h and for nearby dwarf irregulars 
jFerrero et al.ll201 1 ). The evidence, however, is indirect, since halo 
masses are estimated by extrapolating data that probe only the inner 
few kiloparsecs, where most baryons reside. 

We explore here the possibility of using satellite galaxies 
to help constrain the virial masses of dwarf galaxies. The or- 
bital motions of satellite companions have often been used to es- 
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timate halo masses (see, e.g., Zaritskv et al.|l 19971: Erickson et al. 


1999 


; McKay et alj 


20021; Prada et al.1 120031; Ivan den Bosch et al. 


2004 


; Brainerd 2005 


:IConrov et al.l2007l:lwoitak & Mamorj2012l). 



but this work has largely been restricted to systems similar to or 
brighter than the Milky Way. This is partly due to the difficulties 
in obtaining redshifts for faint objects. In addition, satellite com- 
panions are less common around dwarf galaxies than around larger 
systems: the number of satellites brighter than a certain fraction 
of the primary luminosity, JV(> L sat JL pT1 ) , declines strongly to - 
ward fainter primaries (e.g. iGuo et aljioTTI i lWang & Whitel20i2h . 
Dwarf galaxy associations do e xist, but only a handful have been 
observed (e.g.. lTulIv etal]|2006h . 

In ACDM, where satellite galaxies are thought to trace the 
substructure of cold dark matter halos, satellite systems are ex- 
pected around all central galaxies, regardless of luminosity. The 
number of satellites, and their dependence on primary mass, should 
just reflect the abundance of substructure modulated by the depen- 
dence of galaxy formation efficiency on halo mass. 

Substructure abundance has been studied extensively through 
numerical simulations, and shown to be nearly invariant with halo 
mass when expr essed as a function of satellite mass normalized to 
that of the host dMoore et alll999l; iKravtsov et~aT1l2004 iGao et all 
120041 : IWang et al]|2012h . This result, together with the strict con- 
straints on galaxy formation efficiency mentioned above, imply that 
satellite number counts provide useful information on the halo mass 
of dwarf galaxies. In particular, the near self-similarity of cold dark 
matter halos provides an instructive test: if satellite galaxies trace 
substructure then the abundance of luminous satellites should also 
be scale-free on scales where galaxy mass and halo mass are related 
by a power law. 

We explore these issues here by identifying primary-satellite 
systems in galaxy catalogues constructed from the Sloan Digital 
Sky Survey and by comparing them with predictions from a semi- 
analytic mock galaxy catalogue based on the Millennium Simula- 
tions. The plan for this paper is as follows. Sec.|2]describes briefly 
the observational and model datasets while Sec.[3]presents our main 
results. We summarize our main conclusions in Sec. [4] 



difference |AVi. . s | < 3V200 are then classified as satellites (see 
Sec. l2.2t . We have checked that our results are not sensitive to vari- 
ations by factors of a few in these thresholds, nor to the addition of 
0.2-0.35 dex scatter to the M*-M2oo relation. 

We app ly volume and edge corrections to our sample in the 
same way as lWang & White! d2012t) . Completeness for SDSS spec- 
troscopic data is estimated to be ~ 90% for apparent r -band 
magnitudes brighter than m r — 17.7 (Blan tonet all 12005). This 
limit applies to satellite and primary galaxies independently. How- 
ever, satellites are fainter than their centrals, having an absolute 
magnitude difference AM = M^ at - M r pri . This means that 
the effective volume where satellites of a given AM are (almost) 
complete varies strongly with Af^", or equivalently, with primary 
mass. For instance, satellites about a magnitude fainter than their 



hosts will have MI 



-20.5 and —15.2 for primaries with 



stellar masses in the ranges 10.5 < log(M* / Mq) < 11 and 
7.5 < log(M» /Mq) < 8, respectively. In our sample, the limiting 
distance where satellites can be found with AM = 1 is roughly 
320 kpc for primaries with 10.5 < log(M*/M Q ) < 11. This 
halves for primaries with M* ~ 1O 1O M0, and drops to 35 Mpc for 
centrals with 7.5 < log(M»/M0) < 8. Despite the smaller vol- 
ume surveyed for low mass galaxies, there are still enough galaxies 
to probe the satellite population of even the faintest primaries we 
considered. 

One might think that the detection of faint objects would 
be aided by selecting satellites from the photometric catalogue; 
which is complete down to ~ 4 magnitudes fainter than the spec- 
troscopic sample. This requires the stacking of primary galax- 
ies and a statistical background subtraction in order to identify 
the excess count corresponding to satellites jLorrimer et al.1 19941; 
i Lares et alj201 lLlLlu et al 1201 lMNierenberg et al 1201 ll : lGuo et all 
boilUStrigari & Wechsledkod ; IWang & WhitefeOll) . However, 
the signal-to-noise to carry out this subtraction is too low for our 
faint primaries (M* rl < 1O 9 M0). Our sample is then selected 
purely from the spectroscopic catalogue. We discuss briefly pos- 
sible biases affecting faint, low surface brightness companions in 
Sec[3] 



2 DATA AND CATALOGUES 

2.1 Satellite and primary galaxies in SDSS/DR7 

We select primary galaxies spanning a wide range of stellar mass, 
7.5 ^ log(M*/Mo) 11, from the spectroscopic New York Uni- 
versity Value Added Galaxy Catalog (NYU-VAGC). This catalogue 
was built on the basis of the s eventh data release of the Sloan Dig- 
ital S ky Survey (SDSS/DR7: iBlanton et al]|2005l ; lAbazaiian et al.1 
2009). Stellar masses are taken directly from the NYU-VAGC cat- 
alog and are estimated by fitting stellar population synthesis models 
to the fc-corrected galaxy colours. They assume a lChabriej I2003I) 
initial mass function. 

We ensure isolation by imposing two conditions: (i) each pri- 
mary must be the brightest of all objects projected within its virial 
radius with line-of-sight velocities differing by less than three times 
the corresponding virial velocity; and (ii) no primary can be located 
within the virial radius of a more massive system. Virial quanti- 
ties are inferred from the stel lar mass, assumin g the abundance- 
matching M„-M2oo relation o flGuo et al.1 d2010l) and cosmological 
parameters consistent with WMAP1 results S}\ = 0.75, Q. m = 
0.25 and H = 73 km/s dSpergel et al]|2003h . Fainter galaxies 
within a projected distance r p < r2oo and a line-of-sight velocity 



2.2 Satellite and primary galaxies in the semi-analytic model 

W e compare our S DSS results with the semi-analytic catalogue 
of IGuo et al. ]]), based on the dark matter only Millen nium- 
II Simulation (hereafter MS-II, iBovlan-Kolchin et alj|201ll) . The 
model parameters are carefully tuned to reproduce the observed 
abundance of low-redshift galaxies over five orders of magnitude 
in stellar mass and 9 mag in luminosity. The semi-analytic data 
contain full 3D velocity and positional information for all galax- 
ies, and thus enables the evaluation of potential biases that may be 
induced by the limited (projected) data available in observational 
surveys. 

Fig. 1 and 9 of IGuo et al.l d201ll) shows that predictions for 
galaxy mass are reliable in simulated halos resolved with at least 
~ 200-300 dark matter particles. This corresponds roughly to 
M200 ~2x 10 9 M o and M, ~ 1O 6 M in MS-II. We therefore 
consider only galaxies with M* ^ 1O 6 M0 in the analysis below. 

Galaxies in the semi-analytic catalogue inhabit dark matter 
halos an d subhalos ide ntified using the SUBFIND group-finder al- 
gorithm lSpri ngeletall200lh . Primary galaxies are the central ob- 
jects of each halo; all other galaxies within the virial radius are 
considered satellites. The catalogue also includes "orphan" galax- 
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Figure 1. Distributions of projected distance (black) and line-of-sight ve- 
locity difference (red) for satellite galaxies in the semianalytic catalogue, 
binned by primary stellar mass. The minimum and maximum mass of each 
bin is quoted in each panel (in units of log(M« rl /Mq)). The distributions 
peak at around the average virial radius and velocity corresponding to a 
given bin (vertical dotted line). The minimum and maximum virial radius 
and virial velocity of the host haloes contributing to each bin is indicated by 
the horizontal line. Notice that, by definition, satellites are located within 
the virial radius of their hosts. Their relative velocity, however, can exceed 
the virial value (see text for more details). 

ies whose subhalos have been disrupted due to numerical resolution 
effects. The catalogue contains more than 157, 000 primary galax- 
ies in the mass range 7 ^ log(M*/Mo) ^ 12 . 

The semi-analytic data can be "projected" to mimic the same 
satellite identification algorithm used for SDSS data (see, e.g., 
IWang & White! 1201 2| for details). Notice that, because of the dif- 
ferent identification criteria applied, the projected satellite and pri- 
mary samples selected from the mock catalogue differ from the 
3D samples (where we use information about the condition as cen- 
tral/satellite object from the SUBFIND catalogues). This enables us 
to calibrate the parameters of the identification procedure in order 
to minimize the contribution of foreground and background objects 
in our primary/satellite sample. 

Fig[TJshow the distributions of projected distance (r p , shown 
in black) and line-of-sight velocity difference (AVj. . s , shown in 
red) between primaries and "true" satellites, grouped in several bins 
of primary stellar mass. Primaries in each bin span a range of virial 
masses given by the M*-M2oo relation and its scatter. For each 
panel, a horizontal line indicates the minimum and maximum virial 
radius/velocity of the host haloes contributing to that bin. All his- 
tograms peak approximately at the average virial radius and virial 
velocity of host halos in the subsample, indicated by the vertical 
dotted line. 

Although by definition r p < r-™ ax = r2oo, the upper bound 
of the velocity difference is less clear, as the escape velocity typi- 
cally exceeds the virial velocity of a halo substantially in the inner 
regions. The red histograms in Fig. [Tj suggest that the large ma- 
jority of true satellites have line-of-sight velocities that differ from 
their primaries by less than ~ 3V200- These considerations justify 



Figure 2. Cumulative number of satellite galaxies within the vi rial radius 
of prim aries in the semi-analytic model of galaxy formation from lGuo et alj 
1201 ll) . Symbols of different colors indicate the average number of satel- 
lites computed after binning primaries by mass in 0.5 dex-width bins of 
M pr ;. Note that, on average, the number of satellites decreases with de- 
creasing primary mass down to Af pr i ~ 10 10 Mq. Below this mass the 
scaled satellite mass function becomes independent of primary mass. 



the choices r™ ax = r 2 oo and A V,™ ax = 3 V200 made to identify 
satellite/primary systems in the observational sample (Sec. l2.lt . 



3 RESULTS 

As discussed in Sec. Q] the galaxy-halo mass relation is expected 
to leave a clear imprint on the abundance of satellites galaxies as a 
function of primary stellar mass. We explore these ideas in Fig. [2] 
using the semi-analytic catalogue described in Sec. 12.21 

Fig. |2] shows, as a function of satellite-to-primary mass ratio, 
the average number of satellites orbiting primaries of different mass 
in the semi-analytic galaxy catalogue. Primaries are binned in log- 
arithmic Mt, bins of 0.5 dex width; the central mass value is quoted 
in the legend. Satellites are identified in 3D, using the full position 
and velocity information available in the catalogue. 

Fig. |2] shows clearly that the average number of satellites of 
bright primaries increases strongly with M» rI . On average, a pri- 
mary as massive as 10 1,5 Mq is surrounded by roughly 10 satel- 
lites more massive than 0.1 M* n . On the other hand, only ~ 40% 
of primaries as massive as the Milky Way (1O 1O,75 M0) have one 
satellite proportionally as massive. The probability of having a 

companion with mJ at /M* n = 0.1 drops further to ~ 10% for 

M pri = 10 io Mq ^ 

Interestingly, this trend does not hold for lower primary 
masses. The satellite abundance, expressed in terms of m% at /M* 11 , 
becomes independent of primary mass in the dwarf-galaxy regime 
(Mf rl < 1O 1O M ). As discussed in Sec.Q] this reflects the feature- 
less power-law scaling between galaxy and halo masses in these 
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scales and is a prime prediction of the model testable by observa- 
tion. 



In order to take into account how projection effects and the 
presence of interlopers may affect this result, we repeat the anal- 
ysis using only projected positions and line-of-sight velocities, as 
described in Sec. 12,21 This enables us to select primaries and satel- 
lites in identical ways for both model and observational datasets. 

Fig. [3] shows, for the mock (left) and SDSS (right) samples, 
the average number of satellites with r-band magnitude difference 
equal to or smaller than AM = M? l[ - M^ at . This is not only 
the simplest observational analogue of the stellar mass ratio but 
also allows a more straightforward derivation of the completeness 
correction for the catalogue. The use of stellar masses would in- 
troduce extra uncertainties that would have to be compensated by 
more conservative cuts to ensure completeness. Working with mag- 
nitudes allows us to increase the effective statistical power of the 
catalogue. As in Fig. [2] we bin primaries according to their stel- 
lar mass, as indicated in the legends. Error bars correspond to 100 
bootstrap resamplings of the data. 

The left panel in Fig. [3] shows that the projected data behave 
similarly to the 3D sample: the abundance of s atellites at given 
AM increases with Mf rl for bright primaries dGuo et all |20 1 ll ; 
IWang & Whitd[2012T) but becomes independent of mass for dwarf 
primaries (Mf rl < 1O 1O M0). The overall behaviour is in strik- 
ingly good agreement with the SDSS satellite abundances, shown 
on the right. Despite the large error bars in the faint-primary bins 
(an unavoidable consequence of the limited effective volume sur- 
veyed by SDSS) the observed trends in satellite count with pri- 
mary stellar mass closely resemble those in the mock catalogue. 
We interpret this result as providing strong evidence in support of 
the nearly power-law dependence of galaxy mass on halo mass on 
dwarf galaxy scales advocated by semianalytic models of galaxy 
formation in the ACDM scenario. 

One concern regarding this interpretation arises from the com- 
pleteness of SDSS spectroscopic data on dwarf galaxy scales. The 
sample is, on average, 90% complete at our apparent magnitude 
limit of m r = 17.7. We partially correct for biases by weight- 
ing satellites and primaries with their FGOTMAIN value, which 
characterizes the complet eness of SDSS in a given region of the 
sky due to fiber collision (Blan tonetaT]|2005h . However, the com- 
pleteness might worsen if objects of low surface brightness, such 
as many dwarfs, are systematically missed. The results in Fig [3] 
would then represent lower-limits to the true satellite abundance, 
but would still provide useful constraints on the halo mass of their 
hosts. We notice, however, that if our results were strongly affected 
by low-surface brightness biases, the good agreement between ob- 
servations and the semianalytic model would be puzzling. Never- 
theless, this point requires further validation once surveys with im- 
proved surface brightness sensitivity become available. 

The deblending of extended objects into multiple spurious 
"galaxies" could potentially affect our results. As a sanity check, 
we have repeated the analysis of Fig. [3] using the NASA Sloan 
AtlasQ (NSA). This catalogue is based on SDSS-DR8 data but 
uses a different background substraction technique that specif- 
ically improves the identification of galaxies over the origi- 
nal SDSS pipeline, especially for lo w-redshift extended objects 
(Blant on"et alj20lH : lGeha et al.l2012l) . This catalogue only include 

2 http://www.nsatlas.org 



galaxies with redshifts z < 0.055 and is almost complete for 
(model) r-band magnitudes m r < 17.2. Once equivalent cuts are 
applied to our sample we find that the abundance of satellites in 
DR7 and NSA are consistent within their error bars; with a small 
upward offset of DR7 compared to NSA. This shift is negligible for 
Mf rl ^ 1O 1O M0 and less than a factor of ~ 2 for fainter primaries, 
which is comparable to the uncertainties estimated by bootstrap of 
the samples. Importantly, the apparent mass-independence of satel- 
lite counts for primaries less massive than MS 11 — 10 10 A/0 is also 
found in the NSA sample. 

In addition to providing hints about the power-law nature of 
the M* -A/200 relation at the low-mass end, satellite number counts 
can also help constrain its slope. Assuming that satellite and pri- 
mary galaxies follow the same M* oc M^ 00 power-law relation 
and that the unevolved subhalo mass function is self-similar, as pre- 
dicted by CDlvffl then the abundance of faint satellites should scale 
roughly with N(> oc ^l a//3) , where /i» = m» at /MS". Notice 
that this behaviour is only expected over the mass ranges where the 
substructure ma ss function is a pow er law, which typically requires 
fl ^ 0.1 (e.g. iGiocoli et alj|2008h . This imposes an upper limit 
/j,* = O.l' 3 ~ 0.003 on the relative mass of the companions where 
we expect iV(> /x*) oc /il a//8 ' to hold. 

In the semianalytic model, a/ fJ ~ 0.95/2.4 ~ 0.4, which is 
roughly in agreement with the slope of the satellite mass function in 
Fig.[2]measured for low-mass companions m * sat /MS rl ^ 0.003. 
This relation is independent of MS provided the simulation re- 
solves satellites differing by three or more orders of magnitude in 
stellar mass with respect to M» rl . This, in our simulations, hap- 
pens at AIS 11 lO lo A/0. Interestingly, the scatter in stellar mass 
at fixed halo mass, which in th e model is ~ 0.35 dex for the low- 
mass objects dGuo et al]|201ll) . does not seem to impact the slope 
derived from the simple arguments given previously. However, this 
would change if the scatter in the Af*-M2oo relation were strongly 
correlated to the halo mass. In our case this correlation is rather 
mild. 

We conclude from Fig. [3] that the good agreement in shape, 
normalization and slope between SDSS primaries and the mock 
catalogue strongly favours a power-law relation with a steep slope 
Af* oc M200 between stellar mass and halo mass of dwarf galax- 
ies. This agrees with predictions from the semianalytic model of 
IGuo et all d201l|) and from extrapolations of abundanc e -matching 
studie s (e.g., iMoster et alj|20ict iBehroozi et alj|20ld IGuo et all 

Hok]). 



4 DISCUSSION AND CONCLUSIONS 

We study the abundance of satellites as a function of primary stellar 
mass in the Sloan Digital Sky Survey. Using the SDSS/DR7 spec- 
troscopic sample from the NYU-VAGC catalogue we are able to 
extend previous studies to significantly fainter primaries, MS 11 = 
[1O 7 ' 5 -1O 11 ]M0. In agreement with previous work, we find that 
the abundance of satellites exceeding a given satellite-to-primary 
stellar mass ratio, m% /MS", depends strongly on MS 11 for bright 
primaries. On the other hand, the abundance of satellites around 
dwarf primaries, Af» rl < 10 10 Mq, is approximately independent 
of primary stellar mass. 

3 N(> fj) oc fi a for small fi, where fi = m^/^DM ' s ^ e rat '° °f me 
dark matter masses of subhalo at infall time and host at z = 0. Typically 
a ~ -0.75 iGiocoli et al.l2008tlYang et alj201ll) 
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Figure 3. Average number of satellites per primary with a r-band magnitude difference smaller or equal to AM = M^ at — M r ■ Each curve correspond 
to a given primary stellar mass range, as indicated by the labels. The left panel is for mock data from the semi-analytic catalogue. Primary/satellite galaxies 
are identified in projection, as outlined in Sec 12.11 Error bars show uncertainties from 100 bootstrap resamplings. Note the lack of dependence on M* 1 
for M» rl < lO lo Af0. Right: Same but for galaxies in the SDSS/DR7 spectroscopic catalogue. As in the left panel, isolated dwarfs with stellar mass 
7.5 < log(M£ rl /MQ) < 10 seem to populate similar dark matter halos, A/200 ~ 10 1 " — 10 11 Mq according to the simulations. 



These results are in excellent agreement with predictions of 
semi-analytic models within ACDM. These trends arise from the 
mass invariance of substructure in CDM halos and from the vary- 
ing efficiency of galaxy formation as a function of halo mass. On 
dwarf galaxy scales, where the relation between galaxy mass and 
halo mass is well approximated by a steep power law, the invari- 
ance of satellite abundance with primary mass reflects directly the 
scale-free nature of substructure. Around bright galaxies the scaling 
between galaxy mass and halo mass deviates from a simple power 
law, leading to the observed strong increase of satellite abundance 
with increasing primary mass. 

Some caveats on the statements above need to be mentioned. 
The first relates to the definition of halo mass for a galaxy that has 
entered the virial radius of a larger object. The mass in subhalos 
is in gene ral ill-define d, depending not only on the identification 
algorithm (Kneb sTet alj|201l[), but also on time, du e to tidal strip- 
ping (e.g. iTormen et all 1 19981 ; iKlypin et al] 1 19991 ; lHavashi et ail 
2003). The stellar mass, however, is more resilient to tidal effects 
and remains approximately constant afte r accretion onto the hos t 
dWhite& Reel 1 19781 ; ISales et all |2007| ; IPenarrubia et"al] l2008h . 
Thus for satellite galaxies, the virial mass at the moment of accre- 
tion is more closely related to the stellar mass than their present- 
day dark halo mass. Since the abundance of subhaloes according to 
their infall mass -termed the "unevolved" subhalo mass function- 
is also independent of host halo mass when written as a function 
of the relative mas s between satellite and host mJcc /Mom ( e -g- 
Giocoli et "al] |2008l ; lYangetal.ll201 1), the arguments given above 
remain valid. 

The second caveat involves the possible dependence of the 
M*-M2oo relation on redshift and, related to that, whether satel- 
lites and primaries follow the same relation between stellar mass 
and halo mass. Abundance matching models suggest that the link 
between the stellar mass and halo mass evolves weakly with red- 



shift ^^^^^^W^^^^W^^^^^^^Sc Jine 
| 20ld: iLeauthaud et al.||201 ll : lYang et alj|2012l ; iMoster et alj|2012f 
iBehroozi et alj|2012f) . Given tha t surviving satellites are preferen- 
tially accreted at late times (e.g. iGao et alj2004ISales et alj2007l) 
the M,-M2oo relation for satellites and centrals are expected to be 
similar. In particular, self-consistent semi-analytical modelling of 
galaxies shows only sma ll differences between the two (e.g. see 
Fig. 9 of iGuo et alfeoilh . 

Notice that these arguments do not allow for stripping of the 
stars from satellites. Arguably, this could complicate the evolution. 
We note, however, that numerical models for dwarf galaxies sug- 
gest t hat once stellar stripp ing sets in, total disruption soon fol- 
lows jPenarrubia et al.l 2008). We therefore expect partial stripping 
of stars to have minor effects in large statistical samples. On the 
other hand, several studies have suggested that accounting for to- 
tal disruption of satellites is needed to match observations (e.g. 
IWeinmann et alj2006l : lKimm et alj2009h . We address this by com- 
paring observational results with a semi-analytic model that explic- 
itly treats satellite disruption by tidal forces. 



Under the assumption of a ACDM universe, the good agree- 
ment in shape and normalization between satellite counts in SDSS 
and those in the mock catalogue provides support for a steep stellar- 
halo mass relation for dwarfs, consistent with the A#* oc Af| g 
predicted both by semi-analytic models and by extrapolations of 
current abundance-matching analyses. More definitive constraints 
on the slope of the M*-M2oo relation for dwarf galaxies may come 
from a robust determination of the slope of satellite abundances 
around isolated primaries in tandem with studies of the effect of 
scatter in the stellar mass - halo mass relation. Probing increasingly 
fainter companions in observational surveys of the surroundings of 
isolated dwarfs may prove crucial for this goal. 
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